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Abstract.  Wakes  in  a Coulomb  crystal  are  produced  by  “pushing”  with  radiation  pres- 
sure on  a rotating  spheroidal  cloud  of  laser-cooled  ions.  The  wakes  are  stationary 
in  the  lab  frame  and  are  caused  by  the  interference  of  “drum-head”  type  oscillations. 
Velocity  images  of  these  wakes  are  obtained  directly  through  the  dependence  of  the  ion 
fluorescence  on  Doppler  shifts,  and  new  analytical  calculations  accurately  reproduce 
these  experimental  wake  images.  The  technique  demonstrates  a way  to  excite  and 
study  modes,  that  were  not  accessible  with  previous  techniques. 

1 Work  of  the  U.S.  Government;  not  subject  to  U.S.  copyright. 


INTRODUCTION 

Mode  studies  [1-3]  provide  an  excellent  way  to  diagnose,  control,  and  probe 
cold  ion  clouds.  Previously  at  NIST-Boulder,  electric  fields  (due  to  either  applied 
voltages  or  trap  imperfections)  acting  on  the  entire  cloud  were  used  to  excite  modes 
of  relatively  long  wavelength  on  the  order  of  the  cloud’s  size  [2,3].  In  this  article,  we 
present  results  on  a novel  way  to  locally  excite  a spectrum  of  waves  with  relatively 
short  wavelength  by  “pushing”  with  radiation  pressure  on  a crystalhzed  ion  cloud. 
In  the  experiments  shown  here,  the  waves  interfere  “downstream”  from  the  push  to 
produce  a wake  pattern  that  is  stationary  in  the  lab  frame,  analogous  to  the  wake 
behind  a ship  moving  in  deep  water  [4].  The  technique  offers  a way  to  locally  probe 
and  diagnose  ion  clouds  and  demonstrates  a method  for  studying  waves  that  were 
not  accessible  with  previous  techniques.  In  addition,  wakes  in  Coulomb  crystals 
are  also  of  current  theoretical  interest  [5,6],  primarily  due  to  recent  experiments  in 
two-dimensional  (2D)  Dusty  Plasma  crystals  [7,8]. 

EXPERIMENTAL  SETUP 

A schematic  of  the  experimental  setup  [9]  is  shown  in  Fig.  1(a).  We  use  a 
cylindrical  Penning-Malmberg  trap  to  confine  Ntot  — 15000  to  45000  ^Be"''  ions 
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FIGURE  1.  (a)  Schematic  of  the  cylindrical  Penning-Malmberg  trap  and  imaging  diagnostics. 

(b)  Top-view  image  of  the  differential  fluorescence  A/(r,  6)  for  an  n = 0.042  Be+  ion  crystal.  The 
annular  region  defined  by  the  dotted  circles  is  analyzed  to  obtain  dispersion  relationship  data. 

(c)  Side  view  image  of  an  q = 0.042  crystal. 


at  a density  of  n ~ 2 x 10*cm“*  with  an  axial  magnetic  field  of  B = 4.465  T 
(giving  a cyclotron  frequency  of  Q.J2t(  = 7.6  MHz)  and  an  electrostatic  potential 
of  Vg  = — lOOOV  (giving  an  axial  frequency  of  uj^/2'k  = 800kHz).  The  radial 
extent  of  our  clouds  varies  from  Rg  = 0.3  to  2 mm,  which  is  much  smellier  than 
the  radius  of  the  trap  walls  Rg,  = 2.0  cm.  The  ion  clouds  are  spheroidal  in  shape 
(i.e.  an  ellipse  of  revolution)  described  by  an  aspect  ratio  a = Zgj Rg,  where  2Zg 
is  the  axial  extent  of  the  cloud.  The  aspect  ratio  is  directly  related  to  the  rotation 
frequency  which  we  control  with  a dipole  “rotating-wall”  perturbation  [10]. 

Thus,  we  are  able  to  experimentally  set  the  cloud  shape  a by  using  the  rotating 
wall  to  set  uJr-  Here,  uir/27r  ranges  from  42.5  kHz  to  128  kHz  with  a respective  range 
in  aspect  ratio  from  q = 0.005  (corresponding  to  a 2D  single  plane)  to  q = 1.0 
(corresponding  to  a 3D  spherical- shaped  cloud).  For  the  images  in  Fig.  1,  a = 0.042 
with  Wr/27T  = 45  kHz,  and  in  both  images  the  cloud  extended  slightly  beyond  the 
aperture  limit  of  the  cameras. 

The  ions  are  laser-cooled  to  temperatures  T ^ 5 mK  (giving  coupling  parameters 
of  r ^ 200)  using  a tunable  laser  set  slightly  to  the  red  of  a strong  atomic  transition 
in  ®Be"'‘  ions  [9].  The  fluorescence  due  to  the  axial  cooling  beam  (shown  in  Fig.  1 
pointing  in  the  —z  direction)  provides  our  primary  diagnostic.  The  resolution  of  the 
optical  systems  is  4/zm,  whereas,  the  inter-particle  spacing  is  ~ 10 /zm.  Unlike 
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previous  mode  studies  [3],  we  do  not  strobe  the  camera;  instead,  we  simply  collect 
the  fluorescence  continuously  for  30  s to  120  s to  generate  an  image. 

The  relative  intensity  of  the  fluorescence  is  highly  sensitive  to  an  ion’s  axial 
velocity  since  the  cooling  laser  is  tuned  to  a relatively  steep  part  of  the  transition 
curve.  Ions  that  are  moving  towards  the  axial  cooling  beam  (Uj  > 0)  are  Doppler- 
shifted  closer  to  the  resonance  peak,  hence  they  scatter  more  photons  and  fluoresce 
more  strongly;  conversely,  those  moving  away  from  the  beam  {v^  < 0)  are  Doppler- 
shifted  further  from  resonance  and  fluoresce  less  strongly.  For  random  thermal 
motion  the  effects  average  out.  However,  for  coherent  ion  motion  the  changes  in 
fluorescence  enable  the  identification  and  measurement  of  axial  oscillations  using 
this  so-called  “Doppler  velocimetry”  diagnostic  [3]. 

For  the  present  experiments  we  excite  waves  in  our  ion  crystals  by  pushing  on 
them  with  laser  radiation.  The  “push  beam”  is  derived  from  the  same  near-resonsmt 
laser  used  for  cooling;  however,  it  has  a relatively  narrow  waist  (~  50  (im)  and  is 
offset  from  the  rotation  axis  (center  of  the  cloud)  by  =155  pm  to  450  pm. 

EXPERIMENTAL  RESULTS 

Fig.  1(b)  shows  an  example  of  Doppler  velocimetry  for  a laser-induced  wake. 
Here  we  show  the  change  in  fluorescence  AI{r,G)  = I(r,B)  — lo{r,0),  where  we 
subtract  off  a “background”  top-view  image  /o(r,  6)  taken  without  the  push  beam 
from  a top- view  image  I{r,  9)  taken  with  the  push  beam.  The  large  white  spot  is 
due  to  the  push  beam  located  at  a distance  Rpb  « 320  pm,  but  the  alternating  dark 
and  light  arcs  are  variations  in  fluorescence  intensity  due  to  coherent  ion  motion 
[11].  We  estimate  [12]  that  in  this  case  the  peak  change  in  fluorescence  corresponds 
to  a change  in  velocity  by  Avz  ~ ±lm/s.  We  further  estimate  [13]  that  this  Avz 
corresponds  to  a displacement  of  Az  ~ ±0.3  pm,  which  is  much  less  than  the  inter- 
particle spacing.  In  other  words,  we  are  able  to  detect  oscillations  of  very  small 
amplitude  due  to  the  extreme  sensitivity  available  with  our  Doppler  velocimetry 
diagnostic. 

These  laser-induced  wakes  are  analogous  to  wakes  behind  a ship  moving  in  deep 
water  [4].  Due  to  the  radiation  pressure,  ions  receive  a downward  “kick”  as  they 
rotate  through  the  push  beam,  similar  to  the  kick  that  water  experiences  as  it  passes 
under  a moving  ship.  In  both  situations,  the  push  excites  a large  spectrum  of  waves 
traveling  in  all  possible  directions  at  speeds  described  by  a dispersion  relationship. 
The  observed  wakes  are  stationary  in  the  frame  of  the  source  and  are  due  to  the 
constructive  interference  of  waves  that  satisfy  the  stationary  phase  condition  [4]. 
Along  the  direction  of  motion,  this  condition  is  satisfied  by  transverse  waves  with 
a phase  velocity  u)/k  that  matches  the  relative  velocity  v of  the  source.  For  our 
rotating  ion  crystals,  the  motion  is  along  a circle  of  radius  Rrp  with  v =lJt  Rrp. 

We  obtain  dispersion  relationship  data  io{k)  by  assuming  that  the  mode  pattern 
in  an  annular  region  directly  behind  the  push  is  due  primarily  to  transverse  waves; 
we  fit  a damped  [14]  sinusoid  to  the  average  change  in  intensity  in  this  region  to 
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FIGURE  2.  Measured  dispersion  relationship  for  clouds  of  various  aspect  ratios  compared  to  a 
theoretical  dispersion  relationship  for  drum-head  oscillations  on  an  infinite  planar  slab  of  thickness 
2Z. 


obtain  k — 27t/A;  and  use  the  relationship  w/fc  = Wr  Rj^  to  obtain  w.  For  example, 
for  the  annular  region  defined  by  the  dotted  circles  in  Fig.  1(b),  the  distance  between 
peaks  gives  A » 185  pm,  which  in  turn  gives  u)j2'K  « 490  kHz  using  the  calculated 
velocity  Wr  Rpb  ^ 90  m/s. 

Dispersion  relationship  data  obtained  in  this  m2inner  is  shown  in  Fig.  2,  where 
the  different  symbols  correspond  to  different  aspect  ratios.  Here  the  wave  frequency 
is  scaled  by  the  plasma  frequency  Wp  = [2wr(Dc  - and  the  wavenumber  is 

multiplied  by  half  the  cloud  thickness  Z at  the  radial  position  of  the  push  beam. 
For  a single  aspect  ratio,  different  wake  patterns  are  generated  by  changing  Rpb, 
which  effectively  changes  the  relative  velocity  of  the  push  beam. 


THEORY 


To  obtain  a theoretical  dispersion  curve  with  which  to  compare  to  the  data  we 
consider  fluid- like,  drum- head  oscillations  in  an  infinite  planar  slab  of  thickness  2Z. 
Using  a analysis  similar  to  that  in  References  [2]  and  [5],  we  obtain 


tan 


(1) 


For  a given  value  of  kZ  this  equation  has  an  infinite  number  of  solutions;  however, 
we  consider  only  the  lowest-order  mode,  for  which  our  diagnostic  is  most  sensitive 
(for  higher-order  modes  axial  variations  in  fluorescence  will  tend  to  average  out  in 
our  top- view  images).  The  dispersion  relationship  for  these  modes  is  plotted  as  the 
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FIGURE  3.  Experimental  images  of  wakes  for  clouds  with  (a)a  = 0.042  and  (c)o  = 0.25  shown 
alongside  corresponding  theoretical  wakes  (b)  and  (d)  for  the  respective  experimental  parameters. 
The  color  scale  corresponds  to  variations  in  velocity  due  to  coherent  motion  of  the  ions. 


solid  curve  in  Fig.  2.  This  theory  curve  agrees  very  well  with  the  data,  particularly 
for  disk-like  clouds  with  low  aspect  ratio. 

We  are  also  able  to  theoretically  predict  the  experimental  wake  patterns  as  shown 
by  the  comparison  in  Fig.  3.  Here,  the  experimental  image  in  Fig.  3(a)  is  the  same 
as  that  shown  in  Fig.  1(b)  for  an  a = 0.042  cloud,  and  the  one  in  Fig.  3(c)  is  for 
a somewhat  thicker  cloud  with  a = 0.25  and  air/27r  = 60  kHz.  For  the  theoretical 
calculations  (b)  and  (d),  we  take  into  account  the  finite  spot  size  and  sum  over  the 
excited  waves  in  a rotating  slab.  These  calculations  eure  able  to  capture  the  main 
features  of  the  experimental  images  very  well.  The  patterns  are  dominated  by 
arc-shaped  transverse  wakes;  however,  lateral  wakes  also  add  additional  structure. 

We  do  not  observe  Mach  cones  in  our  ion  crystals,  even  under  high  velocity 
conditions  (n  > lO"*  m/s).  This  is  because  the  Mach  condition  is  satisfied  only 
if  the  disturbance  moves  faster  than  all  of  the  excited  waves.  According  to  the 
dispersion  relationship  shown  in  Eq.  1 and  Fig.  2,  there  is  no  theoretical  limit  to 
the  speed  of  the  waves  that  we  consider  here,  i.e.  w/fc  — oo  as  fe  — 0.  Mach  cones 
are  also  not  observed  in  the  case  of  gravity  waves  in  water  for  the  same  reason; 
however,  they  do  occur  for  sound  waves  in  air  and  in  dusty  plasma  crystals  [7,8] 
due  to  the  finite  limit  on  the  speed  of  sound  in  these  systems. 
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SUMMARY  AND  DISCUSSION 


In  summary,  we  use  radiation  pressure  to  generate  stationary  wakes  in  a rotating 
Coulomb  crystal  of  Be"*"  ions.  These  wakes  are  theoretically  well  understood  as 
resulting  from  the  interference  of  drum-head  oscillations.  The  experiments  demon- 
strate a novel  method  of  exciting  waves  in  spheroidal  ion  crystals  with  potentially 
important  applications.  For  example,  by  using  the  dispersion  relationship,  the 
method  can  be  turned  into  a top-view  diagnostic  for  the  cloud  parameters  (i.e.  Ur, 
Ntot,  a,  etc.).  Also,  with  a setup  slightly  different  from  that  presented  here,  it  may 
be  possible  to  study  (as  yet  unobserved)  E x B shear  modes  by  modulating  the 
power  of  a perpendicular  push  beam.  Since  these  modes  occur  only  in  solid-like 
plasmas,  such  experiments  would  provide  a new  tool  with  which  to  study  the  cor- 
relations in  a strongly  coupled  plasma.  In  addition,  these  studies  should  also  be 
useful  in  efforts  to  gain  precise  control  over  the  rotational  motion  of  the  ions,  which 
will  be  necessary  for  planned  experiments  on  quantum  information  processing  using 
ion  crystals  confined  in  a Penning  trap. 

This  research  is  supported  by  the  Office  of  Naval  Research  and  the  National 
Science  Foundation  (Dubin).  We  gratefully  acknowledge  useful  comments  on  this 
manuscript  by  Dr.  Wayne  Itano  and  Dr.  Brana  Jelenkovic. 

REFERENCES 

1.  D.H.E.  Dubin,  Phys.  Rev.  Lett.  66,  2076  (1991). 

2.  J.J.  Bollinger,  D.J.  Heinzen,  F.L.  Moore,  W.M.  Itano,  D.J.  Wineland,  and  D.H.E. 
Dubin,  Phys.  Rev.  A 48,  525  (1993). 

3.  T.B.  Mitchell,  J.J.  Bollinger,  X.-P.  Huang,  and  W.M.  Itano,  Optics  Exp.  2,  314 
(1998). 

4.  Lord  Kelvin  (Sir  W.  Thomson),  Proc.  Royal  Soc.  of  London  42,  80  (1887). 

5.  D.H.E.  Dubin,  Phys.  Plasmas  7,  3895  (2000). 

6.  X.  Wang,  A,  Bhattacharjee,  and  S.  Hu,  Phys.  Rev.  Lett.  86,  2569  (2001). 

7.  See  article  “Waves  in  a 2D  Dusty  Plasma  Crystal”  by  J.  Goree  et.  al.  in  these  ssime 
proceedings. 

8.  D.  Samsonov,  J.  Goree,  Z.W.  Ma,  A.  BhattEicharjee,  H.M.  Thomas,  and  G.E.  Morfill, 
Phys.  Rev.  Lett.  83,  3649  (1999). 

9.  For  more  detail  see  article  “Experimented  Dynamics  of  Stressed,  Strongly  Correlated 
Magnetized  Plasm2is”  by  T.B.  MitcheD,  et.  al.  in  these  same  proceedings. 

10.  X.-P.  Huang,  J.J.  Bollinger,  T.B.  Mitchell,  and  W.M.  Itano,  Phys.  Rev.  Lett.  80,  73 
(1998). 

11.  The  thin  concentric  dark  and  light  circles  are  an  artifact  of  the  crystal  structure. 

12.  For  this  estimate  we  assumed  that  Doppler  broadening  was  negligible  and  used  the 
slope  of  the  transition  curve  at  a detuning  of  10  MHz. 

13.  For  this  estimate  we  used  Az  = Av^/u>  with  u)  = (27t)  490  kHz. 

14.  In  these  experiments  the  waves  are  damped  primarily  by  the  axial  cooling  beam. 


126 


